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ABSTRACT 



The origin, reliability and dispersion of the E pi - £ iso and other spectral energy correlations is a highly debated topic in GRB 
astrophysics. GRB080916C, with its huge radiated energy (E iso ~ 10 55 erg in the 1 keV - 10 GeV cosmological rest-frame energy 
band) and its intense GeV emission measured by Fermi, gives us a unique opportunity to further investigate this issue. We also include 
in our analysis another extremely energetic event, GRB 090323, more recently detected and localized by FermifLAT and showing a 
radiated energy comparable to that of GRB080916C in the 1 keV - 10 MeV energy range. Based on Konus/WJND and Fermi 
spectral measurements, we find that both events are fully consistent with the E pi - E lm correlation (updated to 95 GRBs with the data 
available as of April 2009), thus further confirming and extending it and pointing against a possible flattening or increased dispersion 
at very high energies. This also suggests that the physics behind the emission of peculiarly bright and hard GRBs is the same as for 
medium-bright and soft-weak long events (XRFs), which all follow the correlation. In addition, we find that the normalization of the 
correlation obtained by considering these two GRBs and the other long ones for which E p j was measured with high accuracy by the 
Fermi/GBM are fully consistent with those obtained by other instruments (e.g., BeppoSAX , Swift, KonusfWYND), thus indicating 
that the correlation is not affected significantly by "data truncation" due to detector thresholds and limited energy bands. The very 
recent Fermi/GBM accurate estimate of the peak energy of a very bright and hard short GRB with measured redshift, GRB 090510, 
provides further and robust evidence that short GRBs do not follow the E pl - £ iso correlation and that the E p j - E iso plane can be used 
to discriminate and understand the two classes of events. Prompted by the extension of the spectrum of GRB 0809 16C up to several 
GeVs (in the cosmological rest-frame) without any excess or cut-off, we also investigated if the evaluation of E iso in the commonly 
adopted 1 keV - 10 MeV energy band may bias the E p j - E iso correlation and/or contribute to its scatter. By computing E iso from 1 
keV to 10 GeV, the slope of the correlation becomes slightly flatter, while its dispersion does not change significantly. Finally, we find 
that GRB 0809 16C is also consistent with most of the other spectral energy correlations derived from it, with the possible exception 



of the E a 



- fb correlation. 
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1. Introduction 

Despite the huge observational and theoretical advances in the 
last few years, our understanding of the GRB phenomenon is 
still affected by relevant open issues. Among these, the corre- 
lation between the photon energy at which the vF Y spectrum (in 
the cosmological rest-frame) of the prompt emission peaks, Z? Pi i, 
and the total radiated energy computed by assuming isotropic 
emission, Ei so , in long GRBs is one of the most debated and in- 
triguing. Discovered in 200 2 based on a sampl e of BeppoSAX 



GRBs with known redshift (lAmati et al. 2002I) . the E p>i - E 
correlation was then confirmed and shown to hold for all GRBs, 
soft or bright, with known z and constrained values of E p and 
fiuence, with the only exception of the peculia r sub-energetic 
GRB 980425 d Amati 2006al: lAmatietal. 2007b . The existence 
of such a correlation was also supposed by Lloyd et al. (2000) 
based on the analysis of a sample of bright BATSE GRBs with- 
out measured redshift. The implications of this observational 
evidence can include the physics and geometry of the prompt 
emission, the identification and understanding of different sub- 
classes of GRBs (e.g., short, sub-energetic), th e use of G RBs 
for the estimate of cosmological parameters ( Amati 2 006a; 
Ghirlanda et al. 2006t lAmatietal. 2007t lAmatietal. 2008I) . 



Thus, testing the E p j - E\ so correlation and the other "spec- 
tral energy" correlations derived from it, understanding their 
origin and investigating their dispersion and the existence of 
possible outliers is a relevant issue for GRB physics and cos- 
mology. This can be done in three ways: a) by adding new 
data of GRBs with known redshift detected by different instru- 
ments, each one having its own sensitivity and spectral response 
and thus covering different regions of the E v - fiu ence plane 
dAmafi 2006atlAmati 2006& IGhirlanda et al. 20081) : b) by ver- 
ifying its valid ity with large samples of GRBs with no mea - 
sured redshift dGhirlanda et al. 2005al: iGhirlanda et al. 2 008): 
c) by studying th e behaviour in t he E V j - E\ so plane 
of peculiar GRBs dAmati et al. 20071) . Selection effects on 
this correlation have also been i nvestigated with 



con- 



trasting results dBand & Preece 20051: IGhirlanda et aL~2 005a; 
Butl er etal. 20071: IGhirlanda et al. 20081: iButleretal. 2009: 
Shah moradi & Nemiroff 20091) . 

In this article, the £ Pj i - E^o correlation and other 
spectral energy correlations derived from it confront 
the most energeti c GRBs yet dete cted, GRB080916C 
dGreiner et al. 20091: lAbdoetal. 20091) and GRB 090323 
dvanderHorst 20091: iGolenetskii et al. 2009bl) . In particular, 
GRB 0809 16C, with its huge energy release, with the extension 
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Fig. 1. Light curve of the prompt emission of GRB 080916C as 
measured by the Fermi/GBM - n3 detector (-8-1000 keV). The 
horizontal dotted line is the best-fit of the background level as 
measured before and after the GRB. Also shown (vertical dashed 
lines) are the time intervals for which time-resolved spectra 
from 8 keV up to 10 GeV have been reported by Abdo et al. 
(2009). 




Time from trigger [s] 

Fig. 2. X-ray afterglow light curve of GRB 0809 16C as mea- 
sured by the Swift/XRT in 0.3-10 keV. The continuous line 
shows the best-fit double broken power-law; the dashed line 
shows the triple broken power-law obtained by fixing the last 
slope to 2.4 and corresponding to the 90% c.l. lower limit to 
tb (see text). 



of the spectrum of its prompt emission up to tens of GeV without 
any excess or cut-off, and with the accurate measurements of its 
spectral parameters provided by Fermi/GBM and Konus/WIND, 
gives us a unique opportunity of further testing the robustness 
and extension of these correlations and investigating their prop- 
erties. We also update the E p j - Z?j so correlation by including the 
new detected GRBs with known redshift and E p ; u and compare 
the best estimate of its normalization as obtained by using only 
GRBs detected by Fermi/GBM with those estimated with other 
instruments. Our study is based on published spectral results by 
Konus/WIND, Fermi/GBM, Swift and on specific data analysis 
of publicly available data. 



2. Observations and data analysis 

GRB 0809 16C was detected by the Fenni/GBM on 16 
September 2008 at 00:12:45 UT as a long, multi-peak 
struct ured GRB with duration Tgp ~66 s in 50 - 300 
keV dGoldstein & van der Horst 20081) . The light curve of the 
prompt emission as measured by one of the Fermi/GBM Nal 
detectors which triggered the event is shown in Fig. 1. The 
burst was observed also by AGILE (MCAL, SuperAGILE, 
and ACS), RH ESSI, INTEGRAL ( SPI-ACS), Konus/Wmd and 
MESSENGER dHurley et al. 20081) . Remarkably, very high en- 
ergy photons from GRB 0809 16C were detected by Fermi/L AT 
up to ~10 GeV, with more than 145 photons above 100 
MeV and 14 p hotons above 1 GeV (iTajima et al. 2 008: 
lAbdoetal. 20091) . 

Thanks to the prompt dissemination of the Fermi/LAT and 
IPN positions, GRB 0809 16C was followed-up by Swift and 
other ground telescopes, leading to the detection of both the 
X-ray and optical fading counterparts. Of particular interest are 
the afterglow measurements by Swift/XRT and GROND. The X- 
ray afterglow light curve (Fig. 2) shows the canonical shape: a 
steep decay followed by a flat decay and than a steeper power- 
law decay with index ~1.4 and no break up to ~1.3 Ms from 
the GRB onset dStratta et al. 2008b . The optical afterglow light 
curve shows a differe nt behavior, with stro ng evidence of a sim- 
ple power-law decay dGreiner et al. 20091) . 

Later on, a photo metric redshift of 4.3 5±0.15 was reported 
by the GROND team dGreiner et al. 20 09). By combining this 
value with the fluence and spectral parameters of the prompt 
emission provided by Fenni/GBM and Konus/WIND, 0809 16C 
was found to be the most energetic GRB ever, with an E[ so 
of ~4xl0 54 erg in the standard 1-10000 keV cosmological 
rest-frame energy band. Moreover, the joint spectral analysis 
of Fermi/GBM an d LAT data published by the Fermi team 
(Abd oetal. 20 09) showed that the spectrum extends up to ~1- 
10 GeV without any excess or cut-off and that the E lso computed 
by integrating up to 10 GeV is as huge as ~10 55 erg. 

More recently, another very bright GR B 090323 has bee n 
detected and localized by the Fermi/LAT dOhno et al. 2009). 
Follow-up observations were performed by Swift and other 
ground facilities, leadin g to the discovery of the X-ray, optical 
and radio counterpar ts dKennea et al. 2 009: Upd ike et al. 2009c 
lHarrison et al. 20091) . A spec troscopic redshift of 3.57 was 
measured by Gemini south dChornock et al. 20091) . Spectral 
parameters and fluenc e for GRB 090323 were provided by 
both the Fermi/GBM dvan der Horst ~2 009) and Konus/WIND 
dGolenetskii et al. 2009bl) . Based on the spectrum and fluence 
measured by Konus/WIND and the redshift of 3.57 measured by 
Gemini south, it can be found that the E\ so value of this event 
in the 1-10000 keV cosmological rest-frame energy band is 
~4xl 54 erg, thus comparable to that of GRB 0809 16C. No re- 
fined analysis of the VHE emission measured by the LAT from 
this event has been published. 

In our analysis we used results published in the above ref- 
erences and specific data analysis of the public Fermi/GBM and 
Swift/XRT datifl In particular, for GRB 0809 16C we extracted 
the light curve of each Fermi/GBM detection unit which trig- 
gered the event (n3, n4 and bO) by using the gtbin tool included 
in the data reduction and analysis toolfl The Swift/XRT data of 
this burst were processed using the heasoft package (v.6.4). We 

1 The Fermi/GBM data and analysis tools are available at 
ftp://legacy.gsfc.nasa.gov/fermi/, the Swift/XRT data are available at 
http://swift.gsfc.nasa.gov/docs/swift/archive/ 

2 available at http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 
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Table 1. The 25 GRBs with known redshift and measured E p i 
added to the sample of Amati et al. (2008) in our analysis of the 
£ Pi i - Ei so correlation, resulting in a total of 95 GRBs. 



GRB 




[keV] 


[10 52 erg] 


Instrument^ 


Ref.^ 


020127 


1.9 


290+100 


3.5+0.1 


HET 


(1) 


071003 


1.604 


2077±286 


36±4 


KW 


(2) 


080413 


2.433 


584±180 


8.1±2.0 


BAT/WAM 


(3) 


0804 13B 


1.10 


150±30 


2.4±0.3 


BAT 


(4) 


0805 14B 


1.8 


627±65 


17±4 


KW 


(5) 


080603B 


2.69 


376+100 


11 + 1 


KW 


(6) 


080605 


1.6398 


650±55 


24±2 


KW 


(7) 


080607 


3.036 


1691±226 


188+10 


KW 


(8) 


080721 


2.591 


1741±227 


126+22 


KW 


(9) 


080810 


3.35 


1470+180 


45±5 


GBM 


(10) 


080913 


6.695 


710±350 


8.6±2.5 


BAT/KW 


(11) 


080916 


0.689 


184±18 


1.0+0.1 


BAT/GBIVP 


(12) 


081007 


0.5295 


61±15 


0.16±0.03 


GBM 


(13) 


081008 


1.9685 


261±52 


9.5±0.9 


BAT 


(14) 


081028 


3.038 


234±93 


17±2 


BAT 


(15) 


081118 


2.58 


147±14 


4.3+0.9 


BAT/GBM" 


(16) 


081121 


2.512 


871+123 


26+5 


KW 


(17) 


081222 


2.77 


505±34 


30±3 


GBM 


(18) 


090102 


1.547 


1149+166 


22±4 


KW 


(19) 


090328 


0.736 


1028±312 


13+3 


KW 


(20) 


090418 


1.608 


1567±384 


16+4 


BAT/KW 


(21) 


090423 


8.1 


491±200 


11±3 


GBM 


(22) 


090424 


0.544 


273±50 


4.6+0.9 


GBM 


(23) 


0809 16C 


4.35 


2646±566 


380+80 


GBM/KW 


(24) 


090323 


3.57 


1901+343 


410+50 


KW 


(24) 



" Taken from the GRB Table by J. Greiner and references therein 
( http ://w w w. mpe . mpg. de/j cg/grbgen . html ) . 

* Computed in the 1-10000 keV cosmological rest-frame by assum- 
ing a standard ACDM cosmology with H = 70 km s _1 Mpc~', Q.m= 
0.27 and £1 A = 0.73. 

Instrument(s) that provided the spectral parameters and fluence 
used for the computation of E pi and £ iso : HET = HETE-2; BAT 
= Swift/BAT, KW = Konus/WWD; WAM = Suzaku/W AM; GBM = 
Fermi/GBM. 

d References for spectral parameters and fluence: (1) Sakamoto et 
al. (2005); (2) Golenetskii et al. (2007); (3) Ohno et al. (2008); (4) 
Barthelmy et al. (2008a); (5) Golenetskii et al. (2008a); (6) Golenetskii 
et al. (2008b); (7) Golenetskii et al. (2008c); (8) Golenetskii et al. 
(2008d); (9) Golenetskii et al. (2008e); (10) Meegan et al. (2008); (11) 
Palshin et al. (2008); (12) Bissaldi et al. (2008a) and Baumgartner et 
al. (2008); (13) Bissaldi et al. (2008b); (14) Palmer et al. (2008a); 
(15) Barthelmy et al. (2008b); (16) Palmer et al. (2008b) and Bhat et 
al. (2008); (17) Golenetskii et al. (2008g); (18) Bissaldi & McBreen 

(2008) ; (19) Golenetskii et al. (2009a); (20) Golenetskii et al. (2009c); 
(21) Palshin et al. (2009); (22) von Kienlin (2009); (23) Connaughton 

(2009) ; (24) see text. 

' Spectral parameters from Fermi/GBM and fluence from Swift/BAT. 



ran the task xrtpipeline (v.0. 1 1 .6) applying calibration and stan- 
dard filtering and screening criteria^. 

Radiated energies and luminosities are computed by assum- 
ing a standard ACDM cosmology with Ho = 70 km s _1 Mpc -1 , 
£2m= 0.27 and Qa= 0.73. The quoted uncertainties are at 68% 
c.l., unless differently stated. 



see http://swift.gsfc.nasa.gov/docs/swift/analysis/ 



3. The £pi - £i S0 correlation: update and 
comparison among different instruments. 

In Fig. 3 (right panel) we show the E v j - Ei so correlation for long 
GRBs (short GRBs and the peculiar sub-energetic GRB 980425 
are not included) obtained by adding to the sample of 70 events 
of Amati et al. (2008) 25 more GRBs for which measurements 
of the redshift and/or of the spectral parameters have become 
available in the meanwhile (as of April 2009). As in the previ- 
ous evaluations, E[ so is derived in the 1-10000 keV energy band. 
The £p ; and E lso of these events, together with the redshift and 
the relevant references, are reported in Tab. 1 . These values and 
their uncertainties were computed based on published spectral 
parameters and fluences and following the methods and criteria 
reported, e.g., in Amati (2006) and Amati et al. (2008). As can be 
seen, this updated sample of 95 GRBs is fully consistent with the 
£p,i - Zsjso correlation and its dispersion as derived by Amati et al. 
(2008). This is quantitatively confirmed by the fit with both the 
classical x 2 method and with the maximum likelihood method 
adopted, e.g., by Amati (2006) and Amati et al. (2008), which al- 
lows us to quantify the extrinsic scatter of the correlation (<x ext ). 
We obtain a slope m = 0.57±0.01 and a^ 2 of 594, by means 
of a linear fit to the log(£ , p i) vs. log(£'i so ) data points with the^ 2 
method, and m = 0.54+0.03 and cr ext = 0.18+0.02 (68% c.l.) with 
the maximum likelihood method. These values are fully consis- 
tent with those obtained by Amati et al. (2008). 

Butler et al. (2009), in their study of the selection effects, 
claim that the dispersion and significance of the E p j - Ei s0 cor- 
relation in the intrinsic plane is comparable to that of the E p - 
fluence in the observer plane, and that the normalization of the 
Ep i - Eis correlation depends on the instrument used to detect 
GRBs. In Fig. 3 (left panel), we show the distribution of these 
95 GRBs in the E p - fluence observer plane. In order to al- 
low a reliable comparison, the X and Y scales of this plot cover 
the same orders of magnitude as the E p j - E^o plane shown 
in Fig. 3 (right panel). As can be seen, when we move from 
the observer to the intrinsic plane the dispersion of the corre- 
lation between spectral peak photon energy and fluence (radi- 
ated energy) decreases significantly (cr ext from ~0.31 to ~0. 18 
and^- 2 from 3110 to 594), its extension covers more orders of 
magnitudes, and its significance increases (Spearman's p from 
~0.75 to ~0.88). In Fig. 4 we compare the normalization of the 
E D i - Ei S0 correlation obtained with all the most relevant in- 
struments with that derived by Amati et al. (2008). As can be 
seen, no significant (i.e. above ~lcr) change is found. In par- 
ticular, the Ep i - E^o correlation derived using the Swift GRBs 
in which, unlike those considered by Butler et al. (2009), Z? Pi j 
is really measured with BAT (from the official Swift team cata- 
log by Sakamoto et al. 2009 and/or GCNs) or with broad band 
instruments (mainly KonusfWIND), is fully consistent with that 
determined with other instruments. The E p i - Eis correlation 
derived from GRBs detected with Fermi/GBM is fully consis- 
tent with the correlation as determined from other instruments 
with narrower energy ranges. Given the unprecedented broad en- 
ergy coverage of the GBM (from ~8 keV up to more than 30 
MeV), the derived 2? Pj j - Eiso correlation is certainly not affected 
by biases in the estimate of the spectral parameters (the so called 
"data truncation" effect, see. e.g., Lloyd et al. 2000). 

We also tested the effect of the redshift on the E p vs. fluence 
dependence. Starting from the E p vs. fluence data, we derived 
10000 Ep i - Eiso simulated correlations by randomly exchang- 
ing the z values among the 95 GRBs, and we computed for each 
sample the Spearman's correlation coefficient p between the E p j 
and Eiso values so obtained. We found a p distribution fully con- 
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Fig. 3. Location in the E p - fluence (left) and E V j - E- mo (right) planes of the 95 GRBs with firm redshift and E p estimates as of April 
2009 (see text). In both panels the points corresponding to the extremely energetic GRBs 080916C and 090323 are highlighted. In 
addition, in the left panel we mark with red dots those GRBs with spectral parameters and fluence provided by the Fermi/GBM, 
and in the right panel we also show the GRB 080916C point obtained with E lso computed in the 1 keV - 10 GeV cosmological 
rest-frame energy band and the point corresponding to the short GRB 090510. The continuous lines in the right panel correspond 
to the best-fit power-law and the ±2<x dispersion region of the E p j - E^o correlation as derived by Amati et al. 2008. 



sistent with a Gaussian with centroid ~0.75, which is exactly the 
value obtained for the E p vs. fluence correlation in the observer 
plane, with <x~0.035 and extending up to ~0.85. For compari- 
son, the p value of the true E Pj j - E[ so correlation, as we have 
seen, is ~0.88 which is ~3.8 cr from that obtained from the sim- 
ulation and corresponds to a chance probability less than 1 over 
1000, that the true E P j - E\ so correlation is randomly extracted 
from the simulated ones. 

An exhaustive paper devoted to the discussion of the selec- 
tion effects on the E p j - Ei so correlation is in preparation. 

4. Fermi highly energetic GRBs in the £ p> j - E mo 
plane. 

Based on Fermi/GBM, the fluence of GRB 080916C in 8 keV - 
30 MeV was ~1.9xl0~ 4 erg crrT 2 and its time-averaged spec- 
trum in the same ene rgy band can be fit with a Band func- 
tion dBandetal. 19931) with a = -0.91+0.02, /3 = -2.08 +0.06 
and £p = 424+24 keV ( Ivan der Horst & Goldstein""2 008). The 
Konus/WIND team reported, for the 20 keV - 10 MeV energy 
band, a fluence of (1.24+0.17)xl0~ 4 erg cirT 2 , a 256-ms peak 
flux of (1.19+0. 30)xl0~ 5 erg cirT 2 s _1 and a time-averaged 
spect rum with g = - 1 .04+0.06 , $ = -2.26+0.3 and E p = 505+75 
keV dGolenetskii et al. 2008fl . 

By taking into account these fluences and spectral parame- 
ters, with their uncertainties, the redshift, with its uncertainty, 
provided by GROND, and by integrating the cosmological rest- 
frame spec trum in the commonly adopted 1 keV - 10 MeV en- 
ergy band dAmati et al. 2 002; Amati 2006a), we derive the fol- 
lowing values: £ iso = (3.8±0.8)xl0 54 erg, E Pji = 2646+566 keV. 
As can be seen in Fig. 3, with these values the location of 
GRB 0809 16C in the £ p ; - Ei S0 plane is very close to the best-fit 
power-law obtained with the sample of 70 long GRBs consid- 
ered by Amati et al. (2008). This confirms that GRB080916C 
follows the E p j - E\ so correlation and extends its range of va- 
lidity along E iso by a factor of ~2. If GRB 080916C is excluded 



from the fit of the correlation, the values of the parameters and 
their uncertainties do not change significantly with respect to 
those reported in the previous Section, which is the case when 
the softest/weakest events are excluded. This confirms that the 
significance and characterization of the E p j - E m0 correlation do 
not depend on events at the extremes of the ranges of E p and 

The spectral analysis performed by Abdo et al. (2009), shows 
that the time resolved spectra of this event can be fit with the sim- 
ple Band function from ~8 keV up to more than 1 GeV. This im- 
plies that the Ei so above 10 MeV could be not negligible. Indeed, 
by extending the integration up to 10 GeV (cosmological rest- 
frame), and using the f3 value provided by Fermi/GBM (which, 
given the extension of the energy band of this instrument, is ex- 
pected to be more accurate than that provided by Konus/WIND), 
we obtain a value of E[ so of (1.1+0. 2)xl0 55 erg, which is higher 
by a factor of ~2.5. As can be seen in Fig. 3, with this (huge) 
value of E iso , GRB 0809 16C is still consistent with the E Vt \ - 
Ei s0 correlation within 2cr and extends its dynamic range along 
£i so by about half an order of magnitude. 

In Fig. 3 we also show the location in the E p ^ - E[ so plane of 
the other ultra-energetic GRB detected more recently by Fermi, 
GRB 090323. For this event, no refined analysis of the VHE 
emission measured by the LAT has been published, thus no reli- 
able extrapolation and integration of the spectrum up to the GeV 
range can be done. Hence, we restrict the analysis to the stan- 
dard 1 keV - 10 MeV energy band. In addition, the published 
GBM spectral analysis concerns only the first ~70 s of the event 
(which shows a total duration of ~120 s), and thus these data 
do not provide a reliable estimate of E p j and E\ so . By using the 
spectral parameters a = -0.96+"' 2 , /3 = -2.09+° ^, E p = 416+^ 
keV and the fluence of (2.0+0.3) xlQ- 4 erg cm 2 (20 keV - 10 
MeV) provided by Konus/WIND dGolenetskii et al. 2009bl) . to- 
gether with the redshift of 3.57 measured by Gemini south, we 
find £ iso = (4.1±0.5)xl0 54 erg and £ Pii = 1901+343 keV. These 
values are very close to those of GRB 0809 16C and make also 
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Fig. 4. Location in the E p j - Ei so plane of those GRBs with localization and E p j provided by different instruments. The top panels 
show those GRBs whose detection, localization and spectrum were provided by BeppoSAX (left) and HETE-2 (right). The bottom- 
left panel shows those GRBs detected and localized by Swift/BAT and for which E p has been provided either by BAT itself or by 
other instruments (excluding the Fermi/GBM). The bottom-right panel shows those GRBs for which the localization has been 
provided either by Swift or Fermi/LAT and E p [ has been measured by Fermi/GBM (right panel). In all panels, the continuous lines 
correspond to the best-fit power-law and the +2<x dispersion region of the correlation as computed by including all 95 GRBs with 
known z and E v j and the dashed line is the best-fit power-law obtained by considering the plotted points only. 



this event fully consistent with the E p j - E- 1!io correlation. This 
is further evidence that the newly discovered class of extremely 
energetic GRBs follows the correlation. The detection of more 
GRBs with photons at GeV energies (e.g., from Fermi/LAT) will 
strengthen this result. 

The extension of the spectrum of GRB 0809 16C supports the 
possibility that, at least for a fraction of long GRBs, the com- 
monly adopted 1 keV - 10 MeV cosmological rest-frame energy 
band for the computation of E{ so may lead to an underestimate 
of this quantity and be a source of systematics and extra-scatter 
in the E p j - Ei so correlation. To test this, we considered again 
the sample of Amati et al. (2008) plus the 25 GRBs reported in 
Tab. 1. For each event, we re-computed the Ei so value by ex- 
tending the integration up to 10 GeV using the a and B values 
reported in the literature. For those events without a reported 
value of B, e.g. in the case of a fit with a cut-off power-law, 



we adopted a Band function with B = -2.3. The fit with the 
X 2 method provides m = 0.55+0.01 with a best-fit^ 2 of 619, 
while the maximum likelihood method provides m = 0.5 1 +0.03 
and cr ext = 0.18+0.02 (68% c.l.). We conclude that extending the 
computation of E[ so up to 10 GeV slightly flattens the slope of 
the E p \ - E lso correlation but does not significantly change its 
scatter. 

Finally, very recently the Fermi/LAT detected and local- 
ized GeV emission fro m a bright short (~0.5s) GRB 090510 
dOhno & Pelassa 20 09). This event was also det ected by AGILE 
at energies above 100 MeV dLongo et al. 2009[). By c ombining 
the VLT redshift estimate of z=0 903 dRau et al. 20091) and the 
spectral parameters and fluence obtained with the Fermi/GBM 
dGuiriec et al. 2 009). it results that the Ei so and £p i of this event 
are (4+l)xl0 52 erg and 8370+760 keV, respectively. With these 
values GRB 090510 lies in the E p> j - Ei S0 plane significantly 
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above the region populated by long GRBs (Fig. 3, right panel), 
further confirming that short GRBs do not follow the E p j - Ei S0 
correlation. 



5. GRB 08091 6C and other spectral energy 
correlations. 

After the discovery and first studies of the E p j - E^o corre- 
lation, it was pointed out that E p ^ also correlates with other 
GRB intensity indicators, like the peak lu minosity, L p i so , 
dYonetoku et al. 2004 iGhirlanda et al. 2005bb or the ave rage 
luminosity. L^, dLamb et al. 2004t iGhirlanda et al. 20091) . In 
addition, it was found that by including the break time of the 
afterglow light cu rve , tb , either directly (|Liang & Zhang 2 005 ; 
iNava et al. 20061: IGhirlanda et al. 2007b or by us- 
ing it to derive the jet opening angle and thus com- 
pute the collimati o n-corrected rad iated energy E y 
iGhirlanda et al. 2004 iNavaetal. 20061) . the extrinsic scat- 
ter decreases significantly. As discussed, e.g., by Amati 
(2008), given the strong correlation between E 1&0 , L lso and 
£p,iso> the two-parameters spectral energy correlations are 
in fact equivalent. In addition, in the light of the Swift re- 
sults on X-ray afterglow light curves, the measurement of 
tb and its use to der iv e the jet opening angl e are questioned 
(ICampana et al. 2001 IGhirlanda et al. 20071) . It was also 
proposed that the inclusion of the "high signal time scale", 7o.45, 
introduced and used for variability studi es, reduces the dis - 
persion of the E pl - Ei SQ correlation (Firmani et al. 2006), 
but this propert y was not confirmed by later studies 
dRossietal. 20081; ICollazzi & Schaefer 2008b . Finally, there 
is evidence that, at least for a significant fraction of GRBs, 
the correlation between E Vt \ and luminosity also holds for the 
time-resolved spectra of individual events dLiang et al. 2004t 
Firma niet al. 20081 iFrontera et al. in prepj) . 

Given its extreme energetics and the good sampling of its 
optical and X-ray afterglow light curve, GRB 080916C can be 
used to test also these correlations. 

Regarding the E p \-Li SO correlation, from the 256 ms peak 
flux measured by Konus and by assuming the best-fit model of 
the time averaged spectrum we derive L Pi i so = (1.9+0.6)xl0 54 erg 
s _1 . This value, combined with the E p j value of 2646+566 keV 
above derived, gives a data point fully consistent with this corre- 
lation (Fig. 5a). 

Regarding the E p \- L p ,; so - To.45 correlation, from the back- 
ground subtracted 8-1000 keV light curves obtained with the 
two (n3 and n4) Fermi/GBM Nal detectors (see Fig. 1), we es- 
timated r . 45 (19.5+0.6 s for «3 and 19.2+0.5 s for n4) follow- 
ing the same approach followed by Rossi et al. (2008). The re- 
sult is that GRB 0809 16C is consistent also with this correlation 
(Fig. 5a). 

From the time-resolved spectral analysis reported by Abdo 
et al. (2009), an accurate estimate of E p j was obtained. By us- 
ing these results, we computed the Li so for each of the cor- 
responding time intervals and we reconstructed the track of 
GRB080916C in the E p j- L iso plane (Fig. 5b). As can be 
seen, the spectral and luminosity evolution of this GRB is 
fully consistent with the E D j- Lj so correlation, as typically ob- 
served for bright event s dLiang et al. 2004tlFirmani et al. 2 008; 
IFrontera et al. in prep.1) . The slope of the power-law that best fits 
the 6 GRB 0809 16C data points is -0.4, slightly flatter than the 
commonly found value of ~0.5. This is mostly due to the data 
point corresponding to the first time interval (A) (Fig. 1), which 
slightly deviates from the Zip,;- L lso correlation. 



Regarding the E p j - E y correlation, no evidence of a jet 
break tb is found either in the X-ray light curve (see Fig. 2) 
up to the end of the XRT observations (~1.3 Ms from the trig- 
ger) or in the optica l light curve up to the end of the GROND 
dGreiner et al. 2009) observations (~0.5 Ms from the trigger). 
A 90% lower limit of about 0.5 Ms to the jet break time is ob- 
tained from the X-ray light curve when assuming a typical post- 
break slope of 2.4 (see Fig. 2). By adopting a lower limit of 
0.5 Ms to tb , standard assumptions on efficiency of conversion 
of the fireball kinetic energy into radiated energy, on the ISM 
density and pr ofile, and on the ratio between mass los s rate and 
wind velocity (Nav aetal. 20061 ; IGhirlanda et al. 2007b . we ob- 
tain £y>8.8xl0 51 erg in the case of an homogeneous circum- 
burst medium and E y > 1.6xl0 51 erg in the case of a wind 
medium. If tb >1 Ms, these values move to £' r >1.5xl0 52 erg 
(homogeneous circum-burst medium), £' y >2.4xl0 51 erg (wind 
medium). These lower limits take into account both the uncer- 
tainty on Ei w and on the GRB redshift. As can be seen in Fig. 5c, 
the lower limits to E y obtained with tb ~ 1 Ms lie at around 1 .5- 
2cr from the best-fit law reported by Ghirlanda et al. (2007), 
while for tb ~0.5 Ms they are fully consistent with it. 

Intriguingly, we find that GRB 0809 16C is a possible outlier 
to the E p j - E lso - fb correlation (Fig. 5d). Indeed, its deviation 
from the best-fit power-law (as determined by Ghirlanda et al. 
2007) is >~3.5cr for tb > lMs and more than ~2.2cr for tb > 
0.5Ms. 



6. Discussion 

The extreme energetics of GRB 0809 16C and the fact that 
its spectrum follows the simple Band function without any 
break or excess up to several tens of GeVs (in the cosmo- 
logical rest-frame) challenge GRB prompt emission models. 
For instance, Abdo et al. (2009) and Wang et al. (2009) sug- 
gest that the favored emission mechanism is the standard non- 
thermal synchrotron radiation from shock-accelerated electrons 
within a fireball with bulk Lorentz factor Y > ~600 - 1000 
dAbdoetal. 20091; iGreiner et al. 2009tlLi 20091) . Nevertheless, 
the lack of a synchrotron self-Compton component cannot be 
explained by this scenario and Inverse Compton in residual col- 
lisions m aybe needed to explain time delayed GeV photons 
dLi 2009b . The fact that GRB080916C is fully consistent with 
the E p j - Eiso correlation (Sect. 4 and Fig. 3), and with most of 
the correlations derived from it (Fig. 5), further supports the hy- 
pothesis that, despite its huge isotropic-equivalent radiated en- 
ergy and the extension to its emission up to VHE, the physics 
behind the emission of this event is not peculiar with respect to 
less energetic long GRBs and XRFs. We note that the E p j - E iS0 
correlation itself can be explained within the non-thermal syn- 
chrotron radiation scenario, e.g., by assuming that the minimum 
Lorentz factor, y ,„,>,, and the normalization of the power-law dis- 
tribution of the radiating electrons do not vary significantly from 
burst to burst or by imposing limits to the slope of the correlation 
betwee n the fireball bulk Lorentz factor, T, and the burst lumi- 
nosity dLlovdet al. 20001; IZhang & Meszaros 20021) . The con- 
sistency of time-resolved spectra of GRB 080916C with the E„i- 
-luminosity correlation (Fig. 5) confirms that the prompt emis- 
sion is dominated by a single emission mechanism. However, the 
slight deviation from this correlation of the peak energy and lu- 
minosity measured during the first time interval (Fig. 5-b) may 
suggest that during the rise phase of the GRB the main emission 
mechanism is not still fully at work and other mechanisms may 
play a relevant role. 
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Fig. 5. (a) Correlation between £ p ; and L p - lso . The red triangle corresponds to GRB 0809 16C and the red continuous lines correspond 
to the best-fit and 2cr range as determined by Yonetoku et al. (2004). The blue square is the L P; j so of GRB 080916C multiplied by 
T9 f ft and the blue dotted lines correspond to the best-fit and 2<x range of the E p j- L p j so - 7o.45 correlation as determined by Rossi et 
al. (2008). (b) Correlation between E P j and L; so within the GRB (red triangles) based on the time resolved spectral analysis reported 
by Abdo et al. (2009). The letters indicate the corresponding time interval (see Fig. 1). The blue square corresponds to the luminosity 
of the whole GRB as determined from the time-averaged spectrum. The red continuous lines show the best-fit power-law and the 
2cr region of the £p,i(t) - Lj so (t) correlation (from Ghirlanda et al. 2007); the blue dotted line is the best-fit power-law to the 6 
points of GRB 08091 6C. (c) Correlation between E p j and the collimation-corrected radiated energy E 7 . The red continuous lines 
and triangles refer to the homogeneous circum-burst medium case while the blue dotted lines and squares refer to the wind case. 
The lines indicate the best-fit power-laws and 2<x regions as reported by Ghirlanda et al. (2007). For both cases, we report the lower 
limits to E y corresponding to tb > 1 Ms and tb > 0.5 Ms (see text), (d) Correlation between £ Pi i, £i so and tb . The lines correspond 
to the best-fit power-laws and 2cr region as reported by Ghirlanda et al. (2007). The two lower limits correspond to tb > 1 Ms and 
t b > 0.5 Ms (see text). 



In turn, GRB 0809 16C confirms the robustness of the E p j 
- E{ so correlation at least in the range of intrinsically medium- 
bright GRBs and, when integrating the spectrum up to 10 GeV, 
extends it by ~half an order of magnitude along E\ so (Fig. 3). 
The flattening of the slope predicted in some scenarios like, e.g., 
the multiple subjet model by Toma et al. (2005) or an increase 
of the dispersion at very high energies is not observed. 

The above considerations are further supported by the other 
extremely energetic GRB 090323 detected more recently by the 
FermifLAT, which shows E p j and Ei so values similar to those 
of GRB 0809 16C and thus is also consistent with the £ p i - Ei S0 
correlation (Fig. 3). The recent measurement by Fermi of the 



£ Pi i and GeV emission of the short bright GRB 090510, com- 
bined with the redshift measurement by VLT, provides further 
and strong evidence that short GRBs do not follow the corre- 
lation holding for long ones, and that the E P j - E\ so plane is a 
powerful tool to discriminate between the two classes and un- 
derstand their different emission mechanisms and origin. 

As a part of our study, we have shown (Sect. 3) that: z) the 
distribution of the updated sample of 95 long GRBs with firm 
estimates of E p j and z in the E p j - E[ so plane is fully consistent 
with the slope, normalization and dispersion determined based 
on previous samples (Fig. 3); ii) moving from the observer frame 
(E p -fluence) to the intrinsic plane (E p j-E lso ) the dispersion of 
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the correlation decreases and its significance significantly in- 
creases (Fig. 3); Hi) if we randomly re-distribute the redshift 
values among the 95 GRBs of the sample, the E v j vs. Ei so dis- 
tribution is similar to to that of E p vs. fluence; iv) not only all 
Fermi/GSM GRBs but also all the other long GRBs with known 
redshift (except GRB 980425) which have been detected with 
BeppoSAX , HETE-2, and Swift, provide E v \ - E- lso correlations 
that are fully consistent with each other and with the E p j - Ei S0 
correlation as derived by Amati et al. (2008) (Fig. 4). 

All this evidence contrasts the conclusions by Butler et al. 
2009 ) that the E p \ - E\ so correlation is strongly affected by in- 
strumental effects. In addition, the fact that GRBs detected and 
localized in different energy bands and by different instruments 
all follow the E v j - E^o correlation favour spectral energy corre- 
lations not being strongly affected by selection effects introduced 
in the observational process that leads to the redshift estimate. 
An exhaustive analysis of instrumental and selection effects on 
the Ep^ - Eis is under way and will be reported elsewhere. 

The spectrum of GRB 080916C following the Band function 
without any cut-off up to a few tens of GeVs (in the cosmologi- 
cal rest-frame) may suggest that the commonly adopted 1 keV - 
10 MeV energy band is too narrow for a correct computation of 
Eiso, thus biasing the E p \ - E lso correlation. However, our anal- 
ysis reported in Sect. 4 shows that the extension of the energy 
band up to 10 GeV which Ei so is computed has a marginal im- 
pact on the slope and the dispersion of the correlation, further 
supporting its robustness. 

Finally, testing the consistency of very high energy GRBs 
with the £p ; j - E\ so and other spectral energy correlations 
Sect. 5) is important for their potential use for cosmology 
Ghirlanda et al. 20061: lAmatietal. 2008h . Indeed, due to de- 
tectors sensitivity thresholds and possible evolutionary effects, 
more luminous GRBs are those more easily detectable at high 
redshifts (e.g., Amati 2006). Besides the E p j - Ei so corre- 
lation, which is fully satisfied by both GRB 0809 16C and 
GRB 090323, the lack of accurate enough long term monitor- 
ing of the optical afterglow of these events ( IGreiner et al. 2009t 
iKann et al. 20091) prevents a stringent test of the correlations in- 
volving the break time tb . However, we find that GRB 0809 16C 
deviates by more than ~2.5cr from the best-fit of the E p j - £; so 
- t correlation (Fig. 5), suggesting that either the dispersion of 
this correlation is higher than thought before, or it is not satisfied 
at very high energies. This is an important issue, given that, with 
respect to the E p j - E y correlation, the E p i - E lso - % has the 
advantage, like the simple E p j - E[ so correlation, of being model 
independent. 

We expect that, thanks to Fermi and AGILE, the number of 
these extremely bright GRBs selected on the basis of their GeV 
emission will increase in the near future, giving us the possibility 
to better understand the physics of the prompt emission of GRBs 
and to get important clues on the reliability and origin of spectral 
energy correlations. 
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